Global entanglement and quantum criticality in spin chains 
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Entanglement in quantum XY spin chains of arbitrary length is investigated via a recently- 
developed global measure suitable for generic quantum many-body systems. The entanglement 
surface is determined over the phase diagram, and found to exhibit structure richer than expected. 
Near the critical line, the entanglement is peaked (albeit analytically), consistent with the notion 
that entanglement — the non-factorization of wave functions — reflects quantum correlations. Singu- 
larity does, however, accompany the critical line, as revealed by the divergence of the field-derivative 
of the entanglement along the line. The form of this singularity is dictated by the universality class 
controlling the quantum phase transition. 
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Introduction — Quantum entanglement, a term coined by 
Schrodinger, has been recognized over the past decade 
as the central actor in many quantum information pro- 
cessing tasks, such as teleportation |l|, dense coding [2j, 
secure keydistribution [3j, and perhaps quantum com- 
putation 0| . More recently, entanglement has emerged 
as an actor on the nearby stage of quantum many- 
body physics, especially for systems that exhibit quan- 
tum phase transitions [a, IE 13 01 ! where it can play the 
role of a diagnostic of quantum correlations. Quantum 
phase transitions are transitions between qualitatively 
distinct phases of quantum many-body systems, driven 
by quantum fluctuations. In view of the connection be- 
tween entanglement and quantum correlations, one an- 
ticipates that entanglement will furnish a dramatic sig- 
nature of the quantum critical point. From the viewpoint 
of quantum information, the more entangled a state, the 
more resources it is likely to possess. It is thus desirable 
to study and quantify the degree of entanglement near 
quantum phase transitions. By employing entanglement 
to diagnose many-body quantum states one may obtain 
fresh insight into the quantum many-body problem. 

To date, progress in quantifying entanglement has 
taken place primarily in the domain of bipartite sys- 
tems '10 . Much of the previous work on entanglement 
in quantum phase transitions has been based on bipar- 
tite measures, i.e., focus has been on entanglement ei- 
ther between pairs of parties [E |(| or between a part and 
the remainder of a system Q • For multipartite systems, 
however, the complete characterization of entanglement 
requires the consideration of multipartite entanglement, 
for which a consensus measure has not yet emerged. 

Singular and scaling behavior of entanglement near 
quantum critical points was discovered in important work 
by Osterloh and co-workers 0, who invoked Wootters' 
bipartite concurrence [Tlj as a measure of entanglement. 
In the present letter, we apply a recently-developed global 
measure that provides a holistic, rather than bipartite, 
characterization of the entanglement of quantum many- 
body systems. Our focus is on one-dimensional spin 
systems, specifically ones that are exactly solvable and 



exhibit quantum criticality. For these systems we are 
able to determine the entanglement analytically, and to 
observe that it varies in a singular manner near the 
quantum critical line. This supports the view that 
entanglement — the non-factorization of wave functions — 
reflects quantum correlations. Moreover, the boundaries 
between different phases can be detected by the entan- 
glement. 

Global measure of entanglement — To introduce a mea- 
sure for characterizing the global entanglement, con- 
sider a general, rt-partite, normalized pure state: — 
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If the parties are all 
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spin- 1/2 then each can be taken to have the basis {| T), 14 
)}. Our scheme for analyzing the entanglement involves 
considering how well an entangled state can be approxi- 
mated by some unentangled (normalized) state (e.g. the 
state in which every spin points in a definite direction): 
|$) = (g)™ =1 The proximity of I*) to |<&) is captured 

by their overlap; the entanglement of |\&) is revealed by 
the maximal overlap |l2j | 
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c (*) = max | ($|*) | ; 
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the larger A max is, the less entangled is |^). (Note that 
for a product state, A max is unity.) If the entangled state 
consists of two separate entangled pairs of subsystems, 
A max is the product of the maximal overlaps of the two. 
Hence, it makes sense to quantify the entanglement of 
|$) via the following extensive quantity |13| 



S log2 W^-log 2 A max (*) 



(2) 



This normalizes to unity the entanglement of EPR-Bell 
and iV-party GHZ states, as well as gives zero for unen- 
tangled states. Finite- ./V entanglement is interesting in 
the context of quantum information processing. To char- 
acterize the properties of the quantum critical point we 
use the thermodynamic quantity £ defined by 
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£ N = N- l E XoS2 ^>) 



(3) 



where £n is the entanglement density, i.e., the entangle- 
ment per particle. 
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Quantum XY spin chains and entanglement — We con- 
sider the family of models governed by the Hamiltonian 

j=i v 7 

where r measures the anisotropy between x and y cou- 
plings, h is the transverse external field, lying along 
the z-direction, and we impose periodic boundary con- 
ditions. At r — we have the isotropic XY limit (also 
known as the XX model) and at r = 1, the Ising limit. 
All anisotropic XY models (0 < r < 1) belong to the 
same universality class, i.e., the Ising class, whereas the 
isotropic XX model belongs to a different universality 
class. XY models exhibit three phases (see Fig.0: oscil- 
latory, ferromagnetic and paramagnetic. In contrast to 
the paramagnetic phase, the first two are ordered phases, 
with the oscillatory phase being associated with a charac- 
terstic wavevector, reflecting the modulation of the spin 
correlation functions (see, e.g., Ref. 14|). We shall see 
that the global entanglement detects the boundaries be- 
tween these phases. 

As is well known 0, IfH Il5| . the energy eigenprob- 
lem for the XY spin chain can be solved by a Jordan- 
Wigner transformation, via which the spins are recast 



as fermions, followed by a Bogoliubov transformation, 
which diagonalizes the quadratic Hamiltonian. Having 
found the eigenstates, A max of Eq. Q and hence the en- 
tanglement can be found. To do this, we parametrize the 
separable states via 

|#)= ® [sm(&/2)|T>i + e < *«coB(&/2)|i) i ], (5) 

i—l 

where | t/l) denote spin states parallel/antiparallel to 
the z-axis. Instead of maximizing the overlap with re- 
spect to the 2N real parameters {£,i,4>i}, for the low- 
est two states it is adequate to appeal to the trans- 
lational symmetry and reality of the wave functions. 

Thus taking £j = £ and <pi = we make the Ansatz: 
_ ■ L v N y 

|$(£)) = e % 2^i=i a j | || ... |) f or searching the maxi- 
mal the overlap A max (^I') [l^. This form shows that this 
separable state can be constructed as a global rotation 
of the ground state at h = oo, viz., the separable state 
1 1| ... |). The energy eigenstates are readily expressed 
in terms of the Jordan- Wigner fermion operators, and 
so too are the Ansatz states By working in this 

fermion basis we are able to evaluate the overlaps be- 
tween the two lowest states and the Ansatz states. With 
l^o) (l^i)) denoting the lowest state in the even (odd) 
fermion-number sector, we arrive at the overlaps 



<*.(r,fc)l*(0> n [cos^)(r,/ l )cos 2 (e/2) + sin^)(r,/ l ) S in 2 ^/2)cot(4 j ) w /2) 
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tan 20« (r, h) = r sin k^] N /(h-cos k^ N ); 
4 0) (0 = 1, /^(0 = ViVsin(£/2)cos(£/2), (TV even); 
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f%\0 ee cos(C/2), = s/Nsin^/2), (N odd); 



where a = 0, 1 and to e [0, N — 1] is an integer. The 
above results are exact for arbitrary N, obtained with 
periodic boundary conditions on spins rather than the so- 
called c-cyclic approximation [15j. Given these overlaps, 
we can readily obtain the entanglement of the ground 
state, the first excited state, and any linear superposi- 
tion, cosal^o) + sinal^i) of the two lowest states, for 
arbitrary (r, h) and N , by maximizing the magnitude of 
the overlap with respect to the single real parameter £. 

The formulas that we have just established contain all 
the results that we explore in the present letter. By ana- 
lyzing the structure of Eq. ©, we find that the global en- 
tanglement does provide information on the phase struc- 
ture and critical properties of the quantum spin chains. 
Two of our key results, as captured in Figs. ^ and [21 are 



that: (i) although the entanglement itself is, generically, 
not maximized at the quantum critical line in the (r, h) 
plane, the field- derivative of the entanglement diverges as 
the critical line h = 1 is approached; and (ii) the entan- 
glement vanishes at the disorder line r 2 + h 2 = 1 , which 
separates the oscillatory and ferromagnetic phases. 

Entanglement and quantum criticality — From Eq. @J it 
follows that the thermodynamic limit of the entangle- 
ment density is given by 



S(r, h) 



- — - max 
In 2 e 



djj, In [cos 0(n,r, h) cos 2 (£/2) 



+ sm6»(^,r,ft,)sin 2 (£/2)cot7r/i] , (8) 



where tan 29([i, r,h) = r sin 2irfJ,/(h — cos 2~k\l). 

Figure [3 shows the thermodynamic limit of the entan- 
glement density £(r, h) and its /i-derivative in the ground 
state, as a function of h for three values of r, i.e., three 
slices through the surface shown in Fig. As the r = 1 
slice shows, in the Ising limit the entanglement density is 
small for both small and large h. It increases with h from 
zero, monotonically, albeit very slowly for small h, then 
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FIG. 1: (Color online) Entanglement density (upper) and 
phase diagram (lower) vs. (r, h) for the XY model with 
N = 10 4 spins, which is essentially in the thermodynamic 
limit. There are three phases: O: ordered oscillatory, for 
r 2 + h 2 < 1 and r =^ 0; F: ordered ferromagnetic, between 
r 2 + h 2 > 1 and h < 1; P: paramagnetic, for h > 1. As is 
apparent, there is a sharp rise in the entanglement across the 
line h — 1, which signifies a quantum phase transition. The 
arc h 2 +r 2 = 1, along which the entanglement density is zero 
(see also Fig. |5J, separates phases O and F. Along r — lies 
the XX model, which belongs to a different universality class 
from the anisotropic XY model. 



swiftly rising to a maximum at h sa 1.13 before decreasing 
monotonically upon further increase of h, asymptotically 
to zero. The entanglement maximum does not occur at 
the quantum critical point. However, the derivative of 
the entanglement with respect to h does diverge at the 
critical point h = 1, as shown in the inset. The slice at 
r = 1/2 (shifted, for clarity, half a unit to the right) shows 
qualitatively similar behavior, except that it is finite (al- 
though small) at h = 0, and starts out by decreasing to a 
shallow minimum of zero at h = \/l — r 2 . By constrast, 
the slice at r = (XX) starts out at ft. = at a max- 
imum value of 1 — 27c/(7rln2) as 0.159. (where 7c is 
the Catalan constant), the globally maximal value of the 
entanglement over the entire (r, h) plane. For larger h it 
falls monotonically until it vanishes at h — 1, remaining 
zero for larger h. 

We find that along the line r 2 + h 2 = 1 the entan- 



FIG. 2: (Color online) Entanglement density and its h- 
derivative (inset) for the ground state of three systems at 
N — oo. Solid blue line: Ising (r = 1) limit; dashed black 
line: anisotropic (r = 1/2) XY model; dash-dotted red line: 
(r = 0) XX model. For the sake of clarity, the XY-case curves 
are shifted to the right by 0.5, indicated by the green arrow. 
For the r = 1/2, at h 2 + r 2 — 1 the entanglement density 
vanishes, which is a general property for the anisotropic XY 
model. Note that whilst the the entanglement itself has a non- 
singular maximum at h w 1.1 (Ising), h w 1.04 (XY r — 1/2), 
h — (XX), respectively, it has a singularity at the quantum 
critical point at h = 1, as revealed by the divergence of its 
derivative. 



glement density vanishes in the thermodynamic limit. 
In fact, this line exactly corresponds to the boundary 
separating the oscillatory and ferromagnetic phases; the 
boundary can be characterized by a set of ground states 
with total entanglement of order unity, and thus of zero 
entanglement density. The entanglement density is also 
able to track the phase boundary (h — 1) between the 
ordered and disordered phases. Associated with the 
quantum fluctuations accompanying the transition, the 
entanglement density shows a drastic variation across 
the boundary and the field-derivative diverges all along 
h = 1. The two boundaries separating the three phases 
coalesce at (r, h) = (0, 1), i.e., the XX critical point. Fig- 
ures n and |21 reveal all these features. 

The singular behavior of the entanglement density (JSJ 
can be analyzed in the vicinity of the quantum critical 
line, and we find the asymptotic behavior (for r ^ 0) 



88 
dh 



1 



2nr In 2 



bx\h- II, for \h-l\< 1. 



(9) 



From the arbitrary- results (0 we analyze the approach 
to the thermodynamic limit, in order to further under- 
stand connections with quantum criticality. We focus on 
the exponent v, which governs the divergence at critical- 
ity of the correlation length: L c ~ \h — l| _ly . To do this, 
we compare the divergence of the slope d£N/dh (i) near 
h = 1 (at N = oo), given above, and (ii) for large N at the 
value of h for which the slope is maximal (viz. /i ma x,Ar), 
i.e., d£N/dh\ hia ^ N w 0.230r _1 In A^ + const., obtained 
by analyzing Eq. 10 for various values of r. Then, not- 
ing that (27rln2)~ 1 w 0.2296 and that the logarithmic 
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scaling hypothesis |l7j identifies v with the ratio of the 
amplitudes of these divergences, 0.2296/0.230 w 1, we 
recover the known result that v = 1. 

Compared with r ^ case, the nature of the divergence 
of d£/dh at r — belongs to a different universality class: 



l0g 2 (7T/2) 



(h 



(10) 



From this divergence, the scaling hypothesis, and the 
assumption that the entanglement density is intensive, 
we can infer the known result that the critical exponent 
v = 1/2 for the XX model. In keeping with the crit- 
ical features of the XY-model phase diagram, for any 
small but nonzero value of the anisotropy, the critical di- 
vergence of the entanglement derivative is governed by 
Ising-type behavior. It is only at the r = point that 
the critical behavior of the entanglement is governed by 
the XX universality class. For small r, XX behavior ul- 
timately crosses over to Ising behavior. 

As is to be expected, at finite N the two lowest states 
|<3>o) an d l^i) featuring in Eq. © do not spontaneously 
break the Z 2 symmetry. However, in the thermodynamic 
limit they are degenerate for h < 1, and linear combina- 
tions are also ground states. The question then arises as 
to whether linear combinations that explicitly break Z 2 
symmetry, i.e., the physically relevent states with finite 
spontaneous magnetization, show the same entanglement 
properties. In fact, we see from Eq. © that, in the ther- 
modynamic limit, overlaps for both l^o) and l^i) are 
identical, up to the prefactors /jy and fjp. These pref- 
actors do not contribute to the entanglement density and 
the entanglement density is therefore the same for both 
l^o) and l^i). It follows that, in the thermodynamic 
limit, the results for the entanglement density are insen- 
sitive to the replacement of a symmetric ground state by 
a broken-symmetry one. 



Concluding remarks — In summary, we have quantified 
the global entanglement of the quantum XY spin chain. 
This model exhibits a rich phase structure, the qualita- 
tive features of which are reflected by this entanglement 
measure. Perhaps the most interesting aspect is the di- 
vergence of the field-derivative of the entanglement as the 
critical line Qi = 1) is crossed. Furthermore, in the ther- 
modynamic limit, the entanglement density vanishes on 
the disorder line (r 2 + h 2 = 1). The structure of the en- 
tanglement surface, as a function of the parameters of the 
model (the magnetic field h and the coupling anistotropy 
r), is surprisingly rich. 

We close by pointing towards a deeper connection be- 
tween the global measure of entanglement and the corre- 
lations among quantum fluctuations. The maximal over- 
lap can be decomposed in terms of correlation func- 
tions: 

IN N 
A L X = ^ + max-^|(r-o ! i>+^(r-CTi(g)r-a ; J ) + ---|, 

where translational invariance is assumed and the 
Carteisan coordinates of r can be taken to be 
(sin£, 0, cos£). The two-point correlations appearing in 
the decomposition are related to a bipartite measure of 
entanglement, namely, the concurrence, which shows sim- 
ilar singular behavior [(|. It would be interesting to es- 
tablish the connection between the global entanglement 
and correlations more precisely, e.g., by identifying which 
correlators are responsible for the singular behavior in the 
entanglement and how they relate to the better known 
critical properties. 
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